Effective Hamiltonian with holomorphic variables 
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The pure-quantum self-consistent harmonic approximation (PQSCHA) permits to study a quan- 
tum system by means of an effective classical Hamiltonian - depending on h and temperature - 
and classical-like expressions for the averages of observables. In this work the PQSCHA is derived 
in terms of the holomorphic variables connected to a set of bosonic operators. The holomorphic 
formulation, based on the path integral for the Weyl symbol of the density matrix, makes it possi- 
ble to approach directly general Hamiltonians given in terms of bosonic creation and annihilation 
operators. 



I. INTRODUCTION 

The notion of Effective Potential in quantum statistical 
mechanics was introduced by Feynman JlJ by means of a 
variational method for the path integral with imaginary 
time. The method was later on improved by Giachetti 
and Tognetti [9 and Feynman and Kleinert || in such 
a way to account completely for the quantum contribu- 
tion of the harmonic part of the interaction. The new 
formulation, obtained by using the path integral in the 
Lagrangian form, has since been successfully applied to 
several condensed matter system, usually providing an 
excellent agreement between theoretical results and ex- 
perimental data Q). The Lagrangian formulation, how- 
ever, presents some difficulties for treating systems for 
which the kinetic energy and the potential energy are 
not separated: spin systems are remarkable instances of 
such a situation. In these cases, indeed, the Feynman- 
Jensen inequality fails to hold and the variational princi- 
ple cannot be directly applied. A generalization that per- 
mits to overcome these problems has been found by using 
the path integral in the Hamiltonian formulation and the 
improved procedure for determining an effective Hamil- 
tonian has been called the pure-quantum self-consistent 
harmonic approximation (PQSCHA) S. Several appli- 



cations to condensed matter have demonstrated its use- 
fulness Nevertheless, this generalization can not be 
applied in a natural way to field-theory models, which 
are rather described in terms of creation and destruction 
operators: in the bosonic case the natural classical- like 
counterpart of these models requires the so called holo- 
morphic variables z* and z, connected to the creation 
and destruction operators. 

In order to calculate Hamiltonian path integrals we 
need a recipe to associate functions in the phase space 
to quantum observables Q| . This implies that some care 
is required to deal with ordering problems: of course, 
neglecting the cases where non-local topological terms 
are present, the ordering procedure must not affect the 
final results. It may however occur that some particu- 
lar ordering turns out to be more suited to a given ap- 
proach or to a specific approximation, in the sense that it 
yields more directly the relevant results. The PQSCHA 
is a semiclassical expansion and therefore the Weyl order- 
ing with the associated Wigner distribution Q appears 
to be the most convenient choice. Weyl ordering, how- 
ever, presents also some drawbacks: for instance it is 
somewhat cumbersome and its physical meaning is not 
immediately evident. Moreover, one can observe that a 
low-temperature field theory is rather formulated using 
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the normal ordering, since the related creation and de- 
struction operators naturally define the vacuum (possibly 
close to the ground state) and it allows for an expansion 
in terms of Feynman diagrams. On the other hand, the 
effective Hamiltonian approach has its main usefulness 
at intermediate and high temperatures, where perturba- 
tive approaches are useless while quanticity is still signif- 
icant, since it fully accounts for the classical nonlinearity 
of the model under study, while the quantum character 
is accounted for at the one-loop level through suitable 
renormalization coefficients; furthermore, it can be ver- 
ified that, whatever the ordering prescription one starts 
with, the final result is such that the effective Hamilto- 
nian, in the limit of low quantum coupling and/or in the 
limit of high temperature, reduces to the Weyl symbol of 
the original Hamiltonian operator. 

In this paper we present the derivation of the effec- 
tive Hamiltonian in the framework of PQSCHA for Bosc 
systems in terms of the holomorphic variables. In the 
next Section we present the path integral formulation for 
calculating the density matrix, we introduce the Weyl 
symbols and their composition properties. We finally es- 
tablish a relation between Weyl and normal symbols that 
turns out to be particularly useful in the explicit calcula- 
tion of the former ones. The third Section shows the 
formulation of the PQSCHA in holomorphic variables 
both in one and in many degrees of freedom. We in- 
troduce the renormalization coefficients that enter the 
expressions of the average values of the dynamical vari- 
ables that, in turn, determine the effective Hamiltonian. 
The additional diagonalization, necessary in the case of 
many degrees of freedom, is also given; the same Section 
presents and uses the low-coupling approximation, and 
its specialization to deal with translation invariant sys- 
tems, for which much more detailed results can be deter- 
mined; eventually, we give an example of the application 
of the general theory to a model system with an interac- 
tion term expressed by a product of occupation numbers, 
thus demonstrating the advantages of the present formu- 
lation. We finally must add that, for the sake of a more 
clear exposition, we have collected in the Appendices the 
most relevant and less immediate calculations needed for 
the development of the subject. 



II. PATH INTEGRAL FOR THE DENSITY 
MATRIX IN COMPLEX VARIABLES 

In order to fix a well defined set of notations, we col- 
lect in this short section some basic properties of the 
Weyl and normal representations of operators in terms 
of complex variables. 

Consider a system of N bosons and let a = 

{a^}^—x n, = {o£}ju=i,...,.!V be the creation and 

destruction operators satisfying the Bose commutation 
relations. We then introduce the corresponding classi- 
cal coordinates (z* z) = {(•£*, .z^)} =1 N . An operator 



atO can be represented either in the normal ordered or 
in the Weyl ordered form, namely 

6 = J d/i(fc*,fc) e^V'k*" N (fc*,fc) 

= J dp(k*,k) e l ^ k + ** & ) 6(fc* k) . (1) 

Here the notation l q is used to denote the transposed 
matrix and, for TV pairs of complex variables q* = 
{q^,}^=i,...,n and q = {^1^=1,...,^, we have introduced 
the integration measure 

d,(q^q) = ]J^^]J d ^^. (2) 
J - J - 2m J - J - 7r 

The normal and Weyl symbols, O n (z* z) and 0(z*,z), 
are then the inverse Fourier transforms of O n (/c* k) and 
0(k*, k), according to the general definition 

f(z*, z)= J dp(k*,k) e< z * k + tk * z ) />* k) . (3) 

It is easily verified that 

0(z*,z) = e-¥ dz * dz O n (z*,z) . (4) 

For practical calculations it is sometimes useful to repre- 
sent the above differential operator as follows: 

e-\ tdz ' dz ^{z\z) = (^(z*+ V *,z- V )} a2=i , (5) 

where (( • )) a 2 is the Gaussian average 

((■}}^=^|w, I 7)(-)e-V^ 2 . (6) 

Weyl symbols have the remarkable property (see Ap- 
pendix A) 

Tr(d!0 2 ) = j dp(z*,z) 1 (z*,z)0 2 (z*,z) , (7) 

that does not extend to products of more than two op- 
erators. In particular, if p = is the density opera- 
tor at the equilibrium temperature T — , we have a 
classical-like expression of the statistical averages: 

<0)-^Tr(/50) = i J dfi(z*,z)p(z*,z)0(z*,z) . 

(8) 

If, for instance, we consider a single harmonic oscilla- 
tor with Hamiltonian Ji = u(d^d + ^) at inverse tem- 
perature P it is straightforward to verify that the Weyl 
symbol for H is Ti.(z*, z) = loz* 'z, while for the density 
matrix, letting / = (3us/2, we first determine the normal 
symbol p K {z*,z) = exp[— / — 2z*z exp(— /) sinh/] and, 
using Eqs. (g) and (||), we finally get 
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p(z* z) 



1 -2z*ztanhf 
e J 



cosh/ 



(9) 



Let us finally consider the standard form for the func- 
tion p(z*, z) in terms of path-integral: 

p(z*, z)= j T>[z*(u),z(u)] e^*( M )' *(")] , (10) 

where T>[z*(u), z(u)] is the functional measure and 



S[z*(u), z(u)] is an Euclidean action, depending on the 
functions z*(u) and z(u) defined for < u < (3. The ex- 
ternal variables z* , z may explicitly enter the expression 
of S through the boundary conditions on the functions 
z*(u), z{u). In this work, however, we find it convenient 
to use an expression of the path-integral in which the 
boundary variables .z*(0), z(0), and z*({3), z{j3), are in- 
tegrated over and the dependence on z* , z, is explicit in 
the form of S, namely 



S[z*(u),z(u)] = j du\^ %*(it);z(u) - V(tt)£(u) - H(z*(u),z(u)) 



If*. 
2 



{ t z*(Q)z{P)~ t z*{P)z{0))-{\z%P)~z*{Q)]z~ t z*[z{P)~z{Q)]) . (11) 



For a complete derivation we refer to the exhaustive 
analysis of phase space path- integrals by Berezin ||, 
while a short but self-contained derivation is presented 
in Appendix A. 



III. PQSCHA: ONE DEGREE OF FREEDOM 

Although the main interest of the method is in the ap- 
plications to problems with many degrees of freedom, it 
is much simpler to understand the method if we first de- 



scribe the case of a single degree of freedom. The idea 
is to decompose the path-integral expression for p(z* z) 
into a first sum over all paths sharing the same average 
point, defined as the functional 



1 f 

(z* z) = - J du (z*(u),z(u)) 



(12) 



and a second sum over average points. In order to do 
this we introduce in the path-integral a resolution of the 
identity that fixes the average point to (z* z). It is then 
natural to define the reduced density, 



p(z*, z; 2* z) = j V[z*(u),z{u)]s(z-^ £ duz(u)\ 6 (z* - i du z* (u)j e S[z*{u),z{u)] ^ (13) 
that collects all contributions coming from paths with the average point (z* z), so that the full density reads 

p{z*,z)= / d/i(z*,z) p(z*,z;z*z) . (14) 



We take p(z*, z; z*, z) as an unnormalizcd probability distribution in the variables (z* z) and define its normalization 
constant as exp(— f3H e g(z*, z)), so that 



p(z*, z- z* z) = e -#tteff(2* z) z . -* 

The thermal average of an observable O can then be written 



(O) = | / dp(z*,z) 



dp(z*,z) 0{z*,z) V(z*, z; z* z) 



(15) 



(16) 



and it is natural to interpret exp[—(37i c g(z*, z)] as a 
classical-like effective density, whereas the probability 
distribution P(z*, z; 1* z) describes the additional fluc- 
tuations around the point (z* z) . In the classical limit it 



can be seen that P(z*, z; z* z) — * S((z*, z) — (z* z)) and 
exp[— /3H e ff(z* z)] tends to the classical Boltzmann fac- 
tor; it follows that the probability V describes the pure- 
quantum fluctuations of the particle, thus providing a 
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separation between classical-like and pure-quantum con- 
tribution to (O). Observe that Eq. (16) is exact and 
provides an ideal starting point for approximations pre- 
serving the full classical nonlinear contribution. 

The explicit evaluation of the reduced density 
p(z*, z; z*, z) will be done in a self-consistent approxima- 
tion replacing TL (z*(u), z(u)) in the action (|ll) with a 
trial Hamiltonian quadratic in the displacements from 
the average point, but depending upon the average points 
themselves. This means that the corresponding trial ac- 
tion results to be non local and therefore cannot be ob- 
tained starting from any quantum operator: we are thus 
dealing with a larger class of actions. The parameters of 
7io are optimized independently for any value of (z*,z). 
Explicitly we take 

H.o{z*, z; z* z) — w(z*, z) + E(z*, z)(z* — z*)(z — z) 

+\[F(z\z){z - z) 2 +c.c] , (17) 

where E(z*,z), F(z*,z) and w(z*,z) are functions to be 
optimized. Linear terms have obviously been neglected 
since they do not contribute to So. It is also evident 
that a quadratic 7io yields a Gaussian probability dis- 
tribution Vq and it turns out that it is centered at the 
average point. We therefore introduce the fluctuation 
variables (£*, £) = (z*—z*,z—z) and we shall denote by 
double brackets the averages over Vq: 

«•))=/ d/*«* ,0 ( • ) W+T, z+d; z*, z) . (18) 



Observe, in particular, that such averages can be easily 
evaluated when the moments ((£*£)), ((£ 2 )} and ((£* 2 )) are 
known. 

As shown in Appendix B, the explicit result for 
Po{z*, z; z*, z) turns out to be 



po{z ,z;z ,z) 



/ r -0w 2 
sinh/ ' £(/) 



-2 



(19) 



Here (£*, £) are connected to (£*, £) by a Bogoliubov trans- 
formation, 



(20) 



that diagonalizes the quadratic term, + \(F£ 2 + 

c.c.) — oj^d . £(/) is the Langevin function 

£(/)=coth/-l//, (21) 

where again / = f3u>/2 while u) is given by 



u) 2 (z*,z) = E 2 (z*,z) - \F(z*,z)\ 



(22) 



Since no ambiguity can anymore arise, in the following 
we shall suppress the bar over z* and z. 

The optimization of the parameters w(z*, z), E(z*, z), 
and F(z*, z) is now in order. According to the PQSCHA 
method we require that Vo gives equal averages for the 
original and the trial Hamiltonian, as well as for their 
second derivatives [Q, namely 



(H(z*+C,z+0) = (H (z*H*,z+0} = w(z*,z) + ±u(z*,z)C{f(z*,z)) , 



(23) 



(d z ,d z H(z*+e,z+0) = (d z ,d z Ho(z*+C, z+0) = E(z*,z) 
(d 2 z H(z*+C;Z+0) = (d 2 Ho(z*H\z+0} = F(z*,z) 
(d 2 z *H(z*+e,z+0) = (d 2 ,H (z*+e,zH)} = F*(z*,z) . 



(24) 



Comparing Eqs. |[|) and (JTsJ) , the PQSCHA effective 
Hamiltonian can be written as 

* * 1 sinh/(z* z) 

n e f[(z ,z) = w(z ,z) + -In — ^ ^ — , (25) 

and eventually Eq. (|l^ ) for thermal averages is approxi- 
mated by 

<0> = \ I dp(z\z) (O(z*H*,z+0) e-WrtW*) . 



(26) 

In terms of the transformed variables the second mo- 
ments of Vo, as defined in Eq. (|l|), are (|*|)) = £(/)/2, 



((£*£*)) — — 0) so that in the original ones their 

expressions turn out to be 

£(/) 



D(z*,z) = ({C0)= E 

C(z\z) = ((e)) = -f* 

C*(z*,z) = ((C 2 )) = -F 



2uj ' 
2uj ' 



(27) 



that we call renormalization coefficients. By means of 
them it is useful to introduce a differential operator A as 
follows. Letting 

A(v,v*) = D(v*,v)d z *d z 

+ ±[C(v*,v)d 2 z +C*{v*,v)d 2 ,]. (28) 
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we define the action of any function F(A) of the operator 
A on a function a(z* z) as 



F(A) a(z*, z) = F(A(v*, v)) a(z*, z) 



(v*,v)=(z»,z) 



(29) 



We can then give a useful representation of the double- 
bracket average ((■)). Indeed, the Gaussian smearing of 
a generic function 0(z* z) can be represented as 



(O(z*H\z+0) = e A 0(z*,z) 



(30) 



For practical calculations, the use of the right-hand-side 
of (|30| ) is particularly convenient when 0{z* , z) is a low- 
degree polynomial. 



The function w(z*, z) is determined by Eq. (23), whose 
last term can be expressed as 

\ u{z\ z)C( f(z*, z)) = (e A A) H(z*, z) , (31) 

so that the effective Hamiltonian can be written in a form 
with a more evident renormalization contribution 



H eS (z*,z)= (1-A) 



„ ., * \ 1 , sinh f(z*. z) 
H(z*,z) + - In I ' 
13 }{z*, z) 



(32) 



It appears therefore that, besides the logarithmic term, 
H e s(z*,z) is given by the Weyl Hamiltonian corrected 



by terms of the order of the square of the renormaliza- 
tion coefficients. Moreover, as we shall show in the next 
Section, the form ( |32| ) is an ideal starting point for a fur- 
ther approximation in order to deal with many degrees 
of freedom. 



IV. PQSCHA: MANY DEGREES OF FREEDOM 



In this section we are concerned with the generaliza- 
tion of the method to a system with N degrees of free- 
dom. The self-consistent solution of the counterpart of 
Eqs. (24)and ( f27j ) for arbitrary values of (z* z) is obvi- 
ously rather difficult: a further simplification is there- 
fore in order. This is called the low-coupling approxi- 



mation (LCA) and consists in expanding the parameters 
— which are the now NxN complex matrices E(z*, z), 
F(z*,z), and F*(z*,z) — so to make the renormaliza- 
tion coefficients, and hence also the Gaussian averages 
((•)), independent of the configuration. Eventually, we 
describe this scheme of approximation and we focus, in 
particular, on translation invariant systems. 

Let us consider therefore a system of N bosons. Again 
the evaluation of p(z*, z; z*, z) is done by replacing H 
with a trial Hamiltonian 7io, 



H {z* z; x* z) = w(z*, z) +' (**-**) E(z*, z) (z-z) + \ [ \z-z) F(z* z)(z-z) + c.c] 



(33) 



quadratic in the displacements from the average point 
(z*,z) = ±jf du(z*(u),z(u)) . (34) 

As already mentioned, E, F, and F* are a Hermitean 
and two symmetric NxN complex matrices, respectively. 

In expressing po we again introduce the pure-quantum 
fluctuation variables (£*, £) = (z* —z*. z—z) and the lin- 
ear canonical transformation 

C = E ( R U(^ *) & + S»k(z*, z) £ k ) 
k 

^ = E ( W> ~ z ) Ck + R»k(z*, z) i k ) (35) 
fe 

to new variables (£ , £) diagonalizes the quadratic term 

E [c E ^ &> + \ F »» £» + c - c - 



"k(z*,z) C k ik ■ (36) 



Here and in the subsequent discussion transformed vari- 
ables are labelled by Latin indices (fc, I, ...); Greek indices 



(/i, v, are used for the original ones. The diagonaliza- 
tion (|36|), that is not possible in the most general case, 
can be performed under suitable constraints on the matri- 
ces .E, i* 1 , and F*, (see e.g. |l(|). However, this problem 
is not extremely relevant in the present context, since 
these matrices are not external data but parameters that 
must be optimized: this constraint just restricts the num- 
ber of independent matrix elements. In addition, in most 
applications the canonical transformation ( |35| ) is mainly 
determined by symmetry considerations. 

In analogy to what we did in the previous section, there 
is no ambiguity in suppressing the bar over (z* z). The 
explicit calculation of po, given in Appendix B, leads to 
the following expression for TC c s- 



H eS (z*,z)=w(z*,z) + ~J2^ 



sinh/ fc (z* z) 
fk(z*,z) 



(37) 



with fk(z*,z) — (3u!k(z*, z)/2, while the moments of 
To are = 8 kl £(/ fe (z* z))/2 and = 

((| fe |j)) = 0. The PQSCHA conditions that determine 
the optimized parameters w(z*,z), E(z*,z), F(z*,z), 
and F*(z*, z), are 
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(H(z*+e, z+O) = (Ho(z*+e, z+O) = w(z* z) + - Y, k Mz*, z) C(f k (z* z)) , 



(38) 



(d^d Zv H(z*+C,z+0)) = (d z *d Zu n (z*+C,z+O)) = E^(z*,z) 
(d Zii d Zu H(z*+e,zH))) = (d z »d Zv H {z*+e,zH)} = F^(z*,z) 
(d z ,d zt H(z*+e,zH)} = (d z *d zt H Q (z*+Z*,zH)} = F^z'z) 



(39) 



Generalizing the procedure of the previous section for the case of one degree of freedom, we define the renormalization 
coefficients 



d^ v {z*,z) = = lj2 k ( R U R »k + s; k s, lk ) c(f k ) 

C IJ , J/ (z*,z) = = -j^ k ( R nkS* vk + R vk S* k ) C(fk) 

c; v (z*,z) = ((e;c)} = \Y, k Ku s »k + KuS„k) c{fk) , 



and, starting from 



A(v*,v) = £ [D^{v\v)d z ,d Zv + l -(C^{v*,v)d Zil d z „ + C; v (v*,v)d z .d zZ ) 



(40) 



(41) 



and defining the action of operator A in analogy with Eq. (29), we find that the Gaussian smearing by double- bracket 
average can be put into the form 



(O(z*+e,z+0) = e A 0(z*,z). 
As in the preceding section, it can be easily seen that the effective Hamiltonian reads 



n eff {z*,z) 



l-A) e A l H(z*,z) +-.E ln 



sinh/ fc (;z*,2:) 



0^ f k (z*,z) 



and the general expression for calculating thermal averages takes the form: 

(o) = \\ &»{z\z) (o(z*+Cz+t)} e -m&{z\z) 



(42) 



(43) 



(44) 



We now describe the low-coupling approximation, 
whose main purpose is to make the averages (( ■ }} in- 
dependent of the configuration jf). The self-consistent 
equations need therefore to be solved only once, with a 
great simplification for implementing the method. One 
can think of several ways to obtain this goal: their choice 
substantially depends on the physics of the problem ||. 
The simplest way, however, consists in expanding the 
matrices E(z*, z), F(z*, z) around a self-consistent min- 
imum (z^Zq) of H cS : 



E(z*,z) = E 
F(z*,z) = F 



6E{z*,z) 
■5F(z*,z) , 



(45) 



where E = E(zq : z ), F = F(zq,z q ), thus using the 
convention to drop the arguments of functions evaluated 
at (zq,z ). As a consequence the frequencies will be ex- 



panded to the first order, uj k (z*, z) = to k + 5uo k (z*,z) 
and for the evaluation of the averages we shall only take 
D and C as renormalization coefficients. It is explic- 
itly shown in Appendix C that the effective Hamiltonian 
becomes 



Wcff (^ 



e A ° H{z* 



z)~\ 



e A oH(z*,z Q ) 



sinh f k 

fk 



(46) 



where f k = /3uj k /2 and 



E 



D„ u d zT d z 



~ (c^d Zfi d z 



c; v d z *o zZ 



(47) 
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The general formula for thermal averages ( j44| ) still holds 
provided that the Gaussian smearing of 0(z*, z) is calcu- 
lated with the LCA renormalization coefficients, i.e. as 
exp(A )O(z*z). 

It often occurs that the indices (i, v, . . . refer to the 
sites of a lattice, whose symmetries can sometimes be 
very helpful in order to simplify the analysis, provided 
that the minimum configuration of TL e g(z*, z) shares the 
same property. The calculations are very easy in the 
case of translation symmetry. In particular, for a one- 
dimensional lattice 



Euv — E 



\i—v i 



F — F 



\i—v ) 



(48) 



where the hermiticity of E and the symmetry of F imply 
that the components E a and F a satisfy 



Performing a Fourier transformation, 

k k 

the left-hand-side of (P6I) becomes 



(49) 



(50) 



E 

fe 

where 



E k Ck^k + \(F k ^-k + c.c.) 



, (51) 



Ek — ^ E a t 



-ika 



-ika 



(52) 



As a consequence of (ft9|) Ek is real and Fk — F_ k . If the 
lattice has more than one dimension the generalization 
is obvious and, in particular, the summation over wave- 
vectors is performed over the first Brillouin zone. If there 
are internal degrees of freedom for each lattice site, the 
Fourier transform must be followed by a diagonalization 
in the internal s pac e. Finally, the diagonalization of the 
quadratic form (|3l]) is completed by a Bogoliubov trans- 
formation 

£,k = AkCk - Bk4~k 

e_ fc - -B-fc*!: + A_ fc *i_ fc , (53) 

where the canonical conditions imply 

A* k A k - B* k B k = 1, A k = A_ k , Bk = B_k, (54) 
and it turns out that 

ui k {A* k A k + B* k B k ) = E k , u k A k B k = F k /2 . (55) 
Finally, 



u k — Ek — F k *F k . 



(56) 



By combining the transformations (|50| )and ( p3[) we ob- 
tain the LCA approximation of the matrices R and S 
defined in @: 



= N-^A* k e ik ^, S; k = -N-^B* k e- ik ^ (57) 



Using then (40) we also can specify the LCA renormal- 



ization coefficients: 



k k 
k 



£(A) 

2uj k 



cos k{jj,—v) . 



(58) 



In many cases only a reduced set of these coefficients will 
be explicitly needed in 7i c ff. For instance, any model 
with on-site nonlinearity entails only on-site coefficients 
D = D^p, C = which are independent of /i due to 
translation invariance. 



V. AN APPLICATION 

As an example of a system that is more conveniently 
treated by means of holomorphic variables, we shall con- 
sider a model Hamiltonian with quartic on-site interac- 
tion, 



E* 

id 



H = - y a: a i+d 



J2{Vin i + V 2 h?) , (59) 



where i runs over sites of a Bravais lattice, d denotes 
the displacements towards the Z nearest-neighbour, hi — 
fijoi, and V\ and V 2 > are constants. One can recognize 
in ( p9| ) the Hamiltonian of the well-known Bose-Hubbard 
model 11 1^| . Its Weyl symbol is easily found to be 



N 

H(z*z)=--V 1 -J2*i*+d 



id 



■E [(^ -V2)z?Z i + V 2 (ztz i f 



(60) 



Ti is invariant under global phase transformations, and 
provided that the LCA effective Hamiltonian shares 
the same property, the translation invariant minimum 
(Zq;, zoi) = (zq, zq) of HcS is degenerate; therefore we can 
choose as representative the real minimum with zo — Zq . 
It is easy to verify that equations (|39| ) for the matrices 
E and F become 

En = [(V 1 -V 2 ) + W 2 (z 2 +D)]6 iS -J2 ^J+d » 

d 

F iS = 2V 2 (z 2 +C)S ii , (61) 

and .Fy = Fy. They yield the following Fourier trans- 
forms: 
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E k = (Vx - V 2 ) + 4V 2 (*o + D )-J2 cos ( k - d ) 

d 

F k = F = 2V 2 {z 2 + C) , 
The renormalization coefficients read now 



(62) 



N 

k 

c = --y 



2w k 



2co k ' 

k 



(63) 



with u>l = El-F 2 . 

Using (f46|), the effective Hamiltonian can be eventually 
written as 



H eS {z*,z) = g{f3)+'H{z*,z) 



(64) 



where 



-V 2 N 



2zl{2D + C) + 2D 2 + C 2 

, 1 V^ln sinll ^ k 



(65) 



is a uniform contribution that affects the partition func- 
tion but cancels out in the expression of thermal averages. 
The (real) translation-invariant stationary points of H. c ff 
are given by 



Z(),l = 



-0,2 



V2-V1 + Z 

2V 2 



2D-C 



(66) 



One can check that when Zq 2 < the minimum is zo.i 
and in the above formulas one should replace zq = Zo,i] 
conversely, for Zq 2 > the minimum is zq, 2 . Since the 
calculation of the renormalization coefficients depends on 
the chosen minimum, both cases can simultaneously oc- 
cur: if this is the case, the relevant minimum to be used 
at a given temperature is the one which gives the smaller 
free energy. 



VI. SUMMARY 

We have presented a construction of the effective 
Hamiltonian suitable for the study of the thermodynam- 
ics of field-theory models, with Hamiltonian expressed in 
terms of creation and destruction operators. The con- 
struction is done by a path-integral approach, in terms 
of the holomorphic variables z* and z. The main result 
we presented consists in Eq. (^) for one degree of free- 
dom, and in Eq. ( p4| ) for many degrees of freedom: in the 
latter case the LCA effective Hamiltonian (^6|) allows one 
to perforrnpractical calculations. As in the previous for- 
mulation |4[] in terms of phase-space conjugate variables, 
p and q, the framework retains all classical nonlinear ef- 
fects embodied in the classical-like formulas for the ther- 
mal averages and the partition function. The quantum 
effects related to the quadratic part of the Hamiltonian 
are completely accounted for, while the pure-quantum 
nonlinearity is treated at one-loop (Hartree-Fock) level. 
The low-coupling approximation permits to actually deal 
with many degrees of freedom, so that our results can be 
used for a direct application to finite-temperature field 
theory. We finally observe that this approach can be 
extended to deal with the thermodynamics of Fermi sys- 
tems [pi. 



APPENDIX A: PATH-INTEGRAL FOR p IN TERMS OF COMPLEX VARIABLES 

We give here a short derivation of the expression for p(z*, z) as a path-integral. Although the subject is rather 
well-known, we have chosen to include it since we want to stress those aspects connected with the Weyl symbols. We 
start from the Trotter's formula for a discrete imaginary time u — * si (i = 1, M, e = M//3) and we take then the 
continuum limit: 

p = e -^= lim (l-sH) M . (Al) 

In order to introduce the Weyl formulation, it is necessary to write the expression of the Weyl symbol of the product 
of M operators. Let us first derive the product rule for two operators O and 0±, i.e. the Weyl symbol 0*Oi(z*, z) 
for OOi . Using Eq. (Q) and the Baker-Campbell-Hausdorff identity, we have 

66 1 = J dfx(k* ,k ) J d/x(fcl,fci) e^Hko+k^ + XK+kDa] e^ko-^lkx)^*^^*^ ^ (A2) 

Using now the expression of the inverse Fourier transform and the representation of the <5-function 

5(k*)d(k) = f dfx(z*,z) e C z * k + tk * z ) , (A3) 



s 



we can determine the product rule for the Weyl symbols that is traditionally called star product: 

o*o x {z*, z) = J d^KM) dKKM) e^fa+fcO + X h l+K>] i ( tklhl - thlkl) 

x J d f i(v* 1 ,v 1 )d f ji(zlz 1 ) e-iCvthx + Xvi) e -*(^fci+%;*i) (v* 1 ,v 1 )0 1 (z* 1 ,z 1 ) . 

2[Xvl-z*)(vi-z 1 )-Xv* 1 -zl){v l -z)] 



= 2 



dfi(v*,vi)dfi(z*,zi) 0{vl,Vx)Ox{z\,Zi) e 



(A4) 



The integration of this expression, using (|A3|), immediately yields Eq. (fjj). By iterating Eq. (A4), we obtain the star 
product of three operators O, 0\, O2'. 



0*0!*0 2 (z*, z) = 2 AN J dfi{vt, Vl ) dfi(zt, Zl ) dfi(v* 2 ,v 2 ) dfi(z*,z 2 ) Ofu^i) O x {zl,z x ) 2 (z* 2 ,z 2 ) 

xe - 2 [K-i'2)(i'i-Zi) - Xv*i-z$)(vi-v 2 )} e -2{%v* 2 -z*)(v 2 -z 2 ) - \vt-zt){v 2 -z)] 

It is now easy to see that in general, if we take O = 1 and therefore O = 1, we can write 

M -1 f AI 



(A5) 



* Oi(z*,z) 

i—1 



/- IVl -1 s 1V1 ~\ 

[] Mvlvi) d^zlzi) Oiizlzi) exp - 2^ [\v*-v* +1 ){v^z t )- \v*-z*)( Vl -v l+1 )]\ , 
-i=i -I i=i J 



(A6) 



where = 2- Letting = z i+d (with Ca/+i = z ) anc ^ Az i = z i — ^i— 1> the previous expression becomes 



* O t (z*,z) 

i—1 



r r M 1 f M 

q2MN J Y[ dfl{C * Xi)d ^ zhZt) o t ( z *, Zl ) exp|2]T[*(C + i+A< + i)C- ^"(Ci+i+Azi+i)] 



(A7) 



The integrations over the variables Ci are performed taking into account a different representation of the 5-function, 
namely 



6(k)6(k*) = I dfi{z*,z) e ±(*z* fe - tfc * z ) 



(A8) 



Since no generality is lost in assuming M to be an even number, we eventually obtain 



A I 
* 

i= 



Jf Oi(z*, z) = 2 MN J f[ dn(z*, Zi ) Oiizlzi) 



x exp < 2 



AI/2 k 



M/2 



fe=l i=l i=l 



(A9) 



Finally, according to Eq. (Al), we set Oj = 1 — e7i ~ e , and in the continuum limit the integrand becomes 

r P/2 pu 

du H(z*(u),z(u)) -2 / / du'[*i:*(2u)i:(2u') - *z* (2u') z{2u)] 



0/2 



+2/ du'[V(/3)z(2u') - t z*{2u l )z{f5)\ +2 / du'[ t i:*(2u / ) 



0/2 



1 - V z(2u')] 



(A10) 



The integration over v! is trivial and, by the change of variable u — ► /3— it, we obtain the path integral given by 
Eqs. @ and @, with the measure T>[z*(u), z(u)] formally defined by 



AI 



■D[z*(u),z(u)} = lim 2 MN Y[d^z*, Zl ) . 



(All) 
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APPENDIX B: THE PQSCHA METHOD 



Let us evaluate the reduced density po (|13|) with the trial Hamiltonian Ho, Eq. (|17|). Representing the S- function 
as in Eq. (AS) one gets 



p (z*,z;z*,z) = / dp(v*,v) 



z v—v z 



(u), z(u)] e 5 °[ z * Z ^ ~ P Io du i z *( u ) v ~ "**(«)] ( ( B1 ) 



where So is given by Eq. ( |ll| ) with Ho- We change the integration variables in the path integral to (£*(«),£(«)) = 
(z*(u)— z* z(u) — z J and we make the Bogoliubov transformation 



£*(u) = R*C(u) + S£(u), £(u) = S*|» + 
with R* R — S* S = 1. This permits to diagonalize the quadratic term 

BrW^) + i[FC 2 (n)+c.c] = w£*(u)£(u) . 



(B2) 
(B3) 



It is straightforward to verify that this transformation preserves at the same time the functional measure and the 
form of the exponent in the path integral. Performing the same transformation on (z* z), [z*, z) and (v*, v) we obtain 

-p {z\z-z\z) = e-P w J d»(v*,v) J 2>[f («),&«)] 

x exp | J du - i [f («) («)] + i [f («) £(u)-£* («) |(u)] - w f (u) £(u) 

-^[r(0)^(/3)-f(/3)^(0)] - [[£*(/3)-r(0)](z-f) ~ (5*-r)[|C9)-|(0)]] | . (B4) 
The shift 

(r («),««)) — (r («) + ^ , ec«) - ^) (B5) 

eliminates the linear term, and the outcoming expression of po contains the path integral of a single harmonic oscillator: 



p (z*,z;z*,z) = e-l 3w f dp(i*,v) e -SV#" / ©[£»,£(«)] 



x exp 



du 



e («)c(u)-r («)$(«) -wr 



f(0)|(/3)-f 08) |(0) 



[cm - f (o)] (5 - 1 + A) - (r - r - 11) [im - |(o)] 



Therefore, using Eq. (g) we obtain 

dp(v* , v) e 



Po(z ,z;z ,z) = e ^ 



■v*v I ' (3uj _ 



cosh / 



->{-»(*•-*•-£)( 



(B6) 



z— z+ 



(3oj 



tanh 



/} > (B7) 



where / = /3u/2. Finally, performing the integration over (v*,v), we arrive at Eq. (|19|). 

Eq. (p3) can be derived from Eq. piT); the latter, in turn, can be easily verified in this way: 



'AH(z*,z) = (AH(z*+e,z+0) = D(z*,z)(d z *d z H(z*H*,z+0) + \[C{z*,z){d 2 z H{z*+C,z+^)) +c.c] . (B8) 



Using Eq. ( |24| ) and Eq. (27) the right hand side becomes 

D(z* z)E(z* z) + | [C(z*, 2)F( 2 * z) + c.c] = [E 2 (z*, z) - F*(z*, z)F(z*, z)] 



£(f(z*,z)) 
2u(z*,z) 



(B9) 



and finally, Eq. (31) follows from Eq. (E2J) 
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In order to deal with the case of many degrees of freedom we use the representation (A8) of the S- function for 
implementing the constraint (|34|), so that the expression for po(z*, z; z*, z) becomes 



p (z*,z;z* z) 



dfj,(v*,v) e 



■z*v- *vz 



Jv[z*(u),z(u)] exp{s [**(«),*(«)] -±J*du[^u)v-S}*z(u)]\, 



(BIO) 



where So is the action (|llj) with the trial Hamiltonian fl33| ) . This expression is evaluated by making the linear canonical 
transformation (|3^) which permits to decouple the path integral into a product of one dimensiona l har monic path 
integrals. The same transformation, done over all the pairs of conjugate variables that appear in ( B10| ), yields the 
form po = exp(— f3w) Ilfc(Po)fc for the reduced density. Taking into account (|3^ ) we can use the result for one degree 
of freedom Eq. ([jl]) and we obtain 



po(z*,Z',z*,z) = e: 



^ W IT ^ ~ 



1 Vsinh/ fe C(f k ) 



exp 



C(fk) 



(Bll) 



where f k = ftuj k /2. 



APPENDIX C: LOW-COUPLING APPROXIMATION 



We give here the steps which lead to the LCA expression (|4E|) for H e s, when the LCA is made by expanding the 
renormalization coefficients around the self-consistent minimum (zq,Zo) of Ti e g. 

We start by expanding Eq. ( (43[) to the lowest order in the difference Suik(z*, z) = uik{z*, z) — u> k , recalling the 
convention to drop the arguments of functions evaluated at (zq, zq). The expansion of the logarithmic term is 



1 -sp sinh/ fc (z*,,z) _ 1 -^p sinh f k 



ft 



fk 



Y^^S^z) 



ft^ ln h 



sinh/ fc r ^ 



[e A A H(z*,z)-e A oA o n(z ,z* ) 



where we have taken into account the fact that Eq. (|3l| ) generalizes to 



}^ u *( z > z ) 2 

k 



H(z*,z). 



Letting A = A + 6 A, we have the expansion 



(1 - A)e A l H(z*, z) = \ (1 - SA)e SA 



e A o H(z*, z) - e A A B H(z*, z) = e A °H(z*,z) - e A A H(z*, z) , 



(CI) 



(C2) 



(C3) 



so that, eventually, Eq. 



is obtained. 
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